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Abstract 

Recent advanced Monte Carlo simulations have not found superconductivity and phase sepa- 
ration in the Hubbard model with on-site repulsive electron-electron correlations. We argue that 
microscopic phase separations in cuprate superconductors and colossal magnetoresistance (CMR) 
manganites originate from a strong electron-phonon interaction (EPI) combined with unavoidable 
disorder. Attractive electron correlations, caused by an almost unretarded EPI, are sufficient to 
overcome the direct inter-site Coulomb repulsion in these charge-transfer Mott-Hubbard insulators, 
so that low energy physics is that of small polarons and small bipolarons (real-space electron (hole) 
pairs dressed by phonons). They form clusters localised by disorder below the mobility edge, but 
propagate as the Bloch states above the mobility edge. I identify the Frohlich finite-range EPI 
with optical phonons as the most essential for pairing and phase separation in superconducting 
layered cuprates. The pairing of oxygen holes into heavy bipolarons in the paramagnetic phase 
(current-carrier density collapse (CCDC)) explains also CMR of doped manganites due to mag- 
netic break-up of bipolarons in the ferromagnetic phase. Here I briefly present an explanation of 
high and low-resistance phase coexistence near the ferromagnetic transition as a mixture of pola- 
ronic ferromagnetic and bipolaronic paramagnetic domains due to unavoidable disorder in doped 
manganites. 
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I. INTRODUCTION 



There are still many theories that attempt to explain the phenomenon of high- 
temperature superconductuctivity in cuprates and other related materials. In general, the 
pairing mechanism of carriers could be not only "phononic" as in the BCS theory fl (left- 
hand upper corner in Fig.flT])) or its strong-coupling bipolaronic extension (right-hand 
upper corner in Fig.flT])) , but also "excitonic" |4|, |5(, "plasmonic" 

aa 

, "magnetic" 

"kinetic" [loj], or purely "coolombic" due to a mirror- nested Fermi surface v\ . The BCS 
theory like any mean-field theory is rather universal, so that it describes well the coopera- 
tive quantum phenomenon of superconductivity even with these non-phononic mechanisms, 
if the coupling is weak (left-hand lower corner in Fig.flT])). The main motivation behind these 
concepts is that high superconducting critical temperature, T c , could be achieved by replac- 
ing phonons in the conventional BCS theory by higher frequency bosonic modes, such as 
plasmons, spin waves (pseudomagnons) , or even by the direct Coulomb repulsion combined 
with unconventional pairing symmetries. 



12] , many authors [13|, [l 



Actually, following original proposal by P. W. Anderson 
assumed that the electron-electron interaction in novel superconductors is very strong but 
repulsive and it provides high T c without any phonons (right-hand lower corner in Fig.flT])). 
A motivation for this concept can be found in the earlier work by Kohn and Luttinger 15( , 
who showed that the Cooper pairing of repulsive fermions is possible. However the same 
work clearly showed that T c of repulsive fermions is extremely low, well below the mK 
scale. Nevertheless, the BCS and BCS-like theories (including the Kohn-Luttinger consid- 
eration) heavily rely on the Fermi-liquid model of the normal state, which fails in many 
high-temperature superconductors. If the normal state is not the Fermi-liquid, then there 
is no direct reason to reject the assumption. In fact there is little doubt that strong onsite 
repulsive correlations (Hubbard U) are an essential feature of the cuprates. Indeed all un- 
doped cuprate compounds are insulators with the insulating gap about 2eV or so. But if 
the repulsive correlations are weak, one would expect a metallic behaviour of a half-filled 
(i-band of copper in cuprates, or, at most, a much smaller gap caused by lattice and spin 
distortions (i.e. due to charge and/or spin density waves jla, Q)- It is a strong onsite 
repulsion of ci-electrons in cuprates which results in their parent insulating "Mott" state. 
When onsite correlations are strong and dimensionality is low, there is an alternative to the 
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FIG. 1: A few theories of high-temperature superconductivity. The highest T c is predicted in the 
BCS to bipolaron crossover |3| for the phonon pairing mechanism (upper half of the diagram). 
Lower half of the diagram represents a number of non-phononic mechanisms of pairing. 



usual Fermi-liquid. In Anderson's resonating-valence-bond (RVB) model [12{| the ground 
state supports "topological solitons" (the so-called spinons and holons), such as occur in 
one-dimensional Hubbard model. Theoretically holons could be paired by a superexchange 



18| 



interaction without any additional glue like phonons or spin-waves 

To discriminate one theory with respect to another one has to rely on experimental facts 
and/or on exact theoretical results. Some variational Monte Carlo (VMC) simulations with 



a (projected) BCS-type trial wave function (see, for example 14| and references therein), and 



a number of other analytical and numerical studies appeared to back up superexchange pair- 
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ing. However recent studies by Aimi and Imada [19f|, using an advanced sign-problem-free 
Gaussian-Basis Monte Carlo (GBMC) algorithm have shown that these variational methods, 
as well as other approximations, overestimated the normal state energy and therefore over- 
estimated the condensation energy by at least an order of magnitude, so that the Hubbard 
model does not account for high-temperature superconductivity. The ground state of the 
model is a normal Fermi liquid with no superconductivity and no stripes. This remarkable 
result is in line with earlier numerical studies using the auxiliary-field quantum (AFQMC) 
[20I and constrained-path (CPMC) 21] Monte-Carlo methods, none of which found super- 
conductivity in the Hubbard model. 

On the other hand compelling experimental evidence for a strong EPI has arrived from 
isotope effects [221 ]. high resolution angle resolved photoemission spectroscopies (ARPES) 



231 ] . a number of optical [24|, neutron-scattering 



25. 



261 ]. tunnelling [27| spectroscopies of 



cuprates, and from recent pump-probe experiments [28(. I have suggested that a strong 
long-range Frohlich EPI is the key to both high-temperature superconductivity and colossal 



magnetoresistance 



Our 
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31 



32 



34 



35 



36| of strongly- 



331 ] and some other studies 
coupled polarons and bipolarons have shown that the long-range discrete Frohlich EPI in 
the presence of the strong Coulomb repulsion does not lead to an enormous enhancement of 



the carrier effective mass characteristic of the Hubbard- 
provide a higher (room temperature) superconductivity 



iolstein model (HHM), and could 



II. PHASE SEPARATION IN CUPRATES 



In the strong-coupling regime, where the BCS electron-phonon coupling constant is rel- 
atively large, A > 1, pairing is individual 2], in contrast with the collective Cooper pairing 
Bipolarons survive even in the normal state above their Bose-Einstein condensation 
temperature representing a simplest "cluster" of carriers. While the Frohlich and Coulomb 



interactions alone could not lead to larger clusters like strings or stripes 401 ] . shorter-range 

421 ] , Jahn- Teller 43j, |44| or a strong nonlinear 



interactions as the deformation |41J], Holstein 



451 ] EPIs could favor bound states of more than two carriers. 

Formation of polaronic clusters can be analytically studied in the strong-coupling regime 
in the framework of a generic "Frohlich- Coulomb" model (FCM) {2], Q, 32]. The model 



Hamiltonian explicitly includes a long-range electron-phonon and the Coulomb interactions 
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as well as the kinetic and deformation energies. The implicitly present large Hubbard U 
term prohibits double occupancy and removes the need to distinguish fermionic spins since 
the exchange interaction is negligible compared with the direct Coulomb and the electron- 
phonon interactions in complex oxides. 

Introducing fermionic, c n , and phononic, d ma , operators the FCM Hamiltonian is written 

as 



n^n' 



T(n - n')c^c n / - ^V c (n - n'^CaC^ev 



- ^2 ^agaim - n)(e a • u m _ n )4c n (rf mQ + d ma ) 

a,nm 

+ + (1) 

ma 

where T(n) is the hopping integral in a rigid lattice, e a is the polarization vector of the 
ath vibration coordinate, u m _ n = (m — n)/|m — n| is the unit vector in the direction from 
electron n to ion m, g a (m — n) is the dimensionless EPI function, and V c (n — n') is the 
inter-site Coulomb repulsion. <7 a (m — n) is proportional to the force / m (n) acting between 
the electron on site n and the ion on m. For simplicity we assume that all phonon modes 
are non-dispersive with frequencies u a and include spin in the definition of n (here h = 1). 
Bipolarons on a two-dimensional lattice of ideal octahedra (that can be regarded as a 



simplified model of the copper-oxygen perovskite layer) have been studied in Ref. 32|]. Due 
to poor screening, the hole-ion interaction was taken as purely coulombic, 



g a (m-n) 



m — n 



where a = x,y, z with I y2 accounting for the experimental fact that c-axis 

(^-polarized) phonons couple to in-plane holes stronger than others. The direct hole-hole 
repulsion is 

VJn - n'l 1 ' 



v^ln-n'l 

so that the repulsion between two holes in the nearest neighbour (NN) configuration is V c . 
The nearest neighbour hopping Tjvjv> the next-nearest neighbour (NNN) hopping across 
copper T NNN and the hopping between the pyramids T' NNN have been included, Fig.(j2J). 
The polaron level shift in this model is given by the lattice sum, 



2 ( 1 h 2 \ „ 

E p = 2k x u q > — ^ + — jg = 31.15^0, 

^-^ \ m — n 4 m — n b / 



(2) 



FIG. 2: Four degenerate in-plane bipolaron configurations A, B, C, and D . Some single-polaron 
hoppings are indicated by arrows. (Reproduced from [32]], (c) IoP, 2002.) 

where the factor 2 accounts for two layers of apical sites, and the in-plane lattice constant is 
a = 1 and u a = ujq. For reference, the Cartesian coordinates are n = (n x + 1/2, n y + 1/2, 0), 
m = (m x , m y , h), and integers. The polaron-polaron attraction is 

, IS a 2 ST^ h2 + ( m ~ n ') ' ( m ~ n ) ,n\ 

^— ' m — rr d m — n \ 6 

m 1 11 1 

Performing the lattice summations for the NN, NNN, and NNN' configurations one finds 
V P h = 1.23 E p , 0.80 E p , and 0.82 E p , respectively. As a result, we obtain a net inter-polaron 
interaction as vnn = V c — 1.23 E p , vnnn = tt| — 0.80 E p , v' NNN = % — 0.82 E p , and the 
mass renormalization exponents (see below) as g 2 NN = 0.38(E p /uj), g 2 NNN = 0.60(E p /uj) 
and (g'NNN) 2 = 0.59(E p /lj). At V c > 1.23 E p , no bipolarons are formed and the system 
is a polaronic Fermi liquid. Polarons tunnel in the square lattice with the renormalised 
hopping integrals t = T NN exp(— 0.38_E p /u;) and t' = T NNN exp(— 0.60E p /u) for NN and 
NNN hoppings, respectively. The polaron mass is m* oc 1/ (t + 2t') . 

If V c < 1.23 E p , then intersite NN bipolarons form. The intersite bipolarons tunnel in the 
plane via four resonating (degenerate) configurations A, B, C, and D, as shown in Fig. (J2j) . 
In the first order of the renormalised hopping integral, one should retain only these lowest 
energy configurations and discard all the processes that involve configurations with higher 

6 



energies. These inter-site bipolarons already move in the first order of the single polaron 
hopping. This remarkable property is entirely due to the strong on-site repulsion and long- 
range electron-phonon interactions that leads to a non-trivial connectivity of the lattice. 
This fact combines with a weak renormalization of if yielding a superlight bipolaron with 
the mass m** oc exp(0.60 E p /uS). We recall that in the Holstein model m** oc exp(2E p /u) [2j]. 
Thus the mass of the small Frohlich bipolaron in the perovskite layer scales approximately 
as a cubic root of that of the Holstein bipolaron. 

At even stronger EPI, V c < 1.16E P , NNN bipolarons become stable. More importantly, 
holes can now form 3- and 4-particle clusters. The dominance of the potential energy over 
kinetic in the transformed Hamiltonian enables us to readily investigate these many-polaron 
cases. Three holes placed within one oxygen square have four degenerate states with the 
energy 2(V C - 1.23E P ) + V c /y/2 - 0.80E P . The first-order polaron hopping processes mix 
the states resulting in a ground state linear combination with the energy £3 = 2.71V C — 



3.2QE p — y/At 2 + t' 2 . It is essential that between the squares such triads could move only 
in higher orders of polaron hopping. In the first order, they are immobile. A cluster of 
four holes has only one state within a square of oxygen atoms. Its energy is E4 = A(V C — 
1.23E p ) + 2{V C /V2-0.80E P ) = 5.41 V c - 6.52E p . This cluster, as well as all bigger ones, are 
also immobile in the first order of polaron hopping. Hence a strong EPI combined with the 
Coulomb repulsion could cause clustering of polarons into finite-size mesoscopic textures. 



Importantly QMC studies of mesoscopic textures [44J including lattice deformations and the 
Coulomb repulsion show that pairs (i.e. bipolarons) dominate over phase separation since 
they effectively repel each other [2] . I would like to stress that at distances much larger than 
the lattice constant the polaron-polaron interaction is always repulsive, and the formation 
of infinite clusters, stripes or strings is prohibited (40|. 

III. PHASE SEPARATION IN FERROMAGNETIC MANGANITES 

The conventional double-exchange (DEX) model of the ferromagnetism and colossal mag- 
netoresistance (CMR), proposed half a century ago and generalized more recently to include 
;he electron-phonon interaction, is in conflict with a number of contemporary experiments 



38. 



39(. Among those experiments are site-selective spectroscopes, which show that oxy 



gen p-holes are current carriers, rather than d-electrons in ferromagnetic manganites 46]. 
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Also, some ferromagnetic manganites manifest an insulating-like optical conductivity at all 
temperatures, contradicting the DEX notion that their ferromagnetic phase is metallic 4t| . 
CMR is observed in manganese pyrochlores 48] where DEX is non-existent. 

On the other hand, the pairing of oxygen holes into heavy bipolarons in the paramag- 
netic phase (current carrier-density collapse (CCDC)) and their magnetic break-up in the 
ferromagnetic phase is compatible with the above and many other observations explaining 



22 



49, 



50 



51 



52] 



CMR, isotope effects, and pseudogaps observed in doped manganites 
Different from other models CCDC predicted the first-order phase transition, now firmly 
established in single CMR crystals [531 ] . CCDC was directly observed in the Hall effect at 



the ferromagnetic transition 

More recently we have proposed an explanation of high and low-resistive phase mixing 
near the ferromagnetic transition observed in tunnelling 55| and some other experiments as 
the mixture of polaronic ferromagnetic domains and bipolaronic paramagnetic domains due 
to unavoidable disorder in doped manganites [391] . Using the fact that the phase transition 



38 



49j one 



in a homogeneous system is of the first order in a wide range of parameters 
can average the magnetisation, er(T), with the Gaussian distribution of random transition 
temperatures Tas caused by disorder around the experimental Curie temperature Tq to 
obtain 



a(T) 



-erfc 



T-T, 



c 



(4) 



2 v r 

where V is the distribution width and erfc(z) = (2/tc 1 ^ 2 ) J z °° dyexp(— y 2 ). The CCDC with 
disorder, Eq. (j3j) fits nicely the experimental magnetizations near the transition with phys- 
ically reasonable T of the order of 10K, depending on doping, Fig. Hence, we believe 
that the random distribution of transition temperatures with the width T across the sam- 
ple caused by the randomness of the bipolaron binding energy is responsible for the phase 
coexistence near the transition as seen in the tunnelling experiments 551 ] . 

Resistivity of inhomogeneous two-phase systems has to be calculated numerically. Nev- 
ertheless, the comprehensive numerical simulations are consistent with a simple analytical 
expression for the resistivity of the binary mixture, 



p = pl P21 



(5) 



which is valid in a wide range of the ratios p\/ P2 56) . Here p\p is the resistivity of each 
phase, respectively, and v is the volume fraction of the second phase. 
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FIG. 3: CCDC model ( lines) describes the experimental magnetisation (upper panel) and resis 



tivity (lower panels) near the ferromagnetic transition in Lao./Cao.aM^-jjTi^Os (symbols 
the phase coexistence caused by disorder is taken into account. No fitting parameters are used in 
Eq.© but the experimental resistivity well below and well above the transition and the experi- 
mental magnetization. (Reproduced from [39|], (c) American Physical Society, 2006.) 



In the framework of CCDC, the resistivity of the paramagnetic phase is px(T) = 
f(T) exp(A/2/c#T) and the resistivity of the ferromagnetic phase is p2{T) = <P{T), where 
f(T) and <p(T) are polynomial functions of temperature depending on the scattering mech- 
anisms, and A is the bipolaron binding energy. Well below the transition 0(T) can be 
parameterized as 4>(T) = p^ + aT 2 , and f(T) = bT well above the transition, where the tem- 
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perature independent parameters po, a, A/2 and b are taken directly from the experiment. 
The volume fraction v of the ferromagnetic phase is simply the relative magnetization in 
our model, v = c(T), also available from the experiment. As a result, Eq. §5$) provides the 
quantitative description of p(T) in the transition region without any fitting parameters by 
using the experimental resistivity far away from the transition and the experimental magne- 
tization, as shown in Fig. ([3]). Studies of the low-field magnetoresistance of Snii^Sr^MnOs 
(x=0.45) which was sintered at different elevated temperatures followed by fast cooling also 
found very good qualitative agreement with CCDC [58] . 



Finally, our concept of polaronic metal in ferromagnetic manganites [38|, |5ll] has been 
substantiated by the angle-resolved photoemission spectroscopy data for the bilayer man- 
ganite Lai^Srx.SM^Oy, where a polaron metallic state below Tq has been clearly observed 



59]. 



IV. CONCLUSIONS 



For although high-temperature superconductivity has not yet been targeted as 'the shame 
and despair of theoretical physics', - a label attributed to conventional superconductivity 
during the first half-century after its discovery - the parlous state of current theoretical 
constructions has led to a current consensus that there is no consensus on the theory of 
high-T c superconductivity Nevertheless impressive amount of experimental data (for ex- 



ample 



22 



23 



24 



25 



26 



27 



281 ]) and accurate numerical simulations [19| have ruled out 



the simple Hubbard model as an explanation of high-temperature superconductivity (lower 
right-hand corner in Fig. ([I])). Our view, which I have briefly presented here in connection 
with the phase separation, is that the extension of the BCS theory towards the strong in- 
teraction between electrons and ion vibrations (upper right corner in Fig.(pQ)) describes the 
phenomenon naturally. The high temperature superconductivity exists in the crossover re- 
gion of the EPI strength from the BCS-like polaronic to bipolaronic superconductivity as was 
predicted by us ]3| before the discovery ]60] , proposed as an explanation of high T r in cuprates 



61 



63. 



64 



65, 



66 



67]- 



621 ]. Fig.(j4]), and explored in greater detail after the discovery 
Bipolarons also explain CMR 38] and, combined with disorder, the phase separation in 
manganites [39j . The observation of the pseudogap and nodal quasiparticles in colossal mag- 



netoresistive manganites 68| which have been considered as a characteristic feature of the 
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FIG. 4: Dependence of the superconducting critical temperature on EPI coupling constant. The 



shading shows the ( 
(Reproduced from 



aiVpolaronic domain. The dotted line corresponds to the BCS-Eliashberg theory. 



62j|, (c) American Physical Society, 1988. 



copper oxide, further substantiates our analogy between high-temperature superconducting 
and CMR oxides. 
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